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Abstract 

The  transient  behavior  of  circuits  like  switching  power  converters  and 
switched  capacitor  filters  are  expensive  to  simulate  because  these  circuits 
are  clocked  at  a  frequency  whose  period  is  orders  of  magnitude  smaller 
than  the  time  interval  of  interest  to  the  designer.  It  is  possible  to  re¬ 
duce  the  simulation  time  without  compromising  accuracy  by  exploiting 
the  property  that  the  behavior  of  such  a  circuits  in  a  given  high  frequency 
clock  cycle  is  similar,  but  not  identical,  to  the  behavior  in  the  preceeding 
and  following  cycles.  In  particular,  the  "envelope”  of  the  high-freqnency 
clock  can  be  followed  by  accurately  computing  the  circuit  behavior  over 
occasional  cycles.  In  this  paper  the  implementation  of  such  an  envelope¬ 
following  method  that  is  particularly  efficient  for  switching  power  and 
filter  circuits  is  described,  and  results  demonstrating  the  method’s  effec¬ 
tiveness  are  presented. 
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1  Introduction 
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In  general,  analog  circuit  designers  rely  heavily  on  circuit  simulation  programs 
like  SPICE  [nagel75]  or  ASTAP  [weeks73]  to  insure  the  correctness  and  the  per¬ 
formance  of  their  designs.  These  programs  simulate  a  circuit  by  first  construct¬ 
ing  a  system  of  differential  equations  that  describes  the  circuit,  and  then  solving 
the  system  numerically  with  a  time  discretization  method  such  as  backward- 
Euler.  When  applied  to  circuits  like  switching  power  converters  or  switched- 
capacitor  filters,  such  classical  circuit  simulation  algorithms  become  extraordi¬ 
narily  computationally  expensive.  This  is  because  switching  power  converters 
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and  switched-capacitor  filters  use  high  frequency  clocks  whose  periods  are  typi¬ 
cally  orders  of  magnitude  smaller  than  the  time  intervals  of  interest  to  a  designer. 
The  nature  of  the  calculations  used  in  a  circuit  simulator  implies  that  to  con¬ 
struct  the  solution  over  the  time  interval  of  interest,  an  accurate  solution  must 
be  computed  for  every  cycle  of  the  high  frequency  clock  in  the  interval,  and  this 
can  involve  hundreds  of  cycles.  The  infeasibility  of  simulating  such  circuits  with 
classical  techniques  has  led  frustrated  designers  to  develop  specialized  analog 
computers  on  which  to  do  their  simulation  [kassakian79]. 

Fast  approximate  techniques  have  been  developed  for  simulating  switching 
power  converters  and  switched  capacitor  filters,  but  the  two  problems  are  ap¬ 
proached  quite  differently.  The  most  common  approach  to  simulating  switched 
capacitor  filters  is  first  to  break  the  circuit  up  into  functional  blocks  such  as  op¬ 
erational  amplifiers  and  switches,  and  replace  each  with  a  simple  macromodel. 
It  is  then  assumed  that  after  each  clock  transition,  every  node  in  the  circuit 
reaches  its  equilibrium  point  before  another  transition  occurs.  This  assump¬ 
tion,  along  with  the  use  of  algebraic  macromodels,  allows  the  filter  to  be  treated 
as  a  low  order  discrete-time  system  whose  solution  for  hunderds  of  clock  cycles 
can  be  computed  quickly  (tsividis79,  deman80j.  The  most  common  techniques 
for  simulating  switching  power  converters  is  to  treat  the  switches  as  ideal,  and 
the  remaining  circuitry  as  linear.  With  this  approximation,  the  solution  over 
hundreds  of  clock  cycles  can  be  computed  rapidly  [hsiao87]. 

Although  programs  based  on  the  above  techniques  have  served  designers  well, 
they  are  based  on  idealizations  of  the  circuits  involved  which  may  eliminate  be¬ 
havior  that  is  important  to  a  designer.  In  this  paper  we  present  an  approach 
for  the  detailed  transient  simulation  of  switching  power  and  filter  circuits  which 
does  not  involve  any  idealization  of  the  behavior,  and  is  much  more  efficient  than 
classical  direct  methods  when  the  clock  period  is  small  compared  to  the  sim¬ 
ulation  interval.  This  method,  referred  to  as  envelope-following  [petzoldfil] 
exploits  the  property  of  such  circuits  that  the  node  voltage  waveforms  over  a 
given  high  frequency  clock  cycle  are  similar,  but  not  exact  duplicates,  of  the 
node  voltages  waveforms  in  proceeding  or  following  cycles.  This  suggests  that 
it  is  possible  to  construct  a  solution  accurate  over  many  high  frequency  clock 
cycles  by  calculating  the  solution  accurately  for  a  few  selected  cycles. 

In  the  next  section,  we  present  the  details  of  an  envelope-following  method 
that  is  effective  for  many  types  of  switching  circuits.  In  Section  3  we  describe 
some  of  the  computations  involved  in  the  method.  Their  implementation  in 
the  program  Niiswit  along  with  results  from  using  Nitawit  to  simulate  several 
switching  power  and  filter  circuits  is  described  in  Section  4. 
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2  The  Envelope-Following  Method 

Most  circuits  can  be  described  by  a  system  of  differential  equations  of  the  form 

lp(*(t),  «(0)  +  /(*(<),  «(<))=  0.  (1) 

where  z(t)  €  5tN ,  the  state,  is  the  vector  of  capacitor  voltages  and  inductor 
currents,  u(t)  €  is  the  vector  of  input  sources,  p(z(t),u(t))  €  &N  is  the 
vector  of  capacitor  charges  and  inductor  fluxes,  and  f(x(t),u(t))  €  $N  is  the 
vector  of  resistive  currents  and  inductor  voltages.  If  the  state  z  is  known  at 
some  time  t0,  it  is  possible  to  solve  (1)  and  compute  the  state  at  some  later 
time  .  In  general,  one  can  write 

z(ti)  =  z(to)  +  ^(z(t0),to,<i)  (2) 

where  *  :  Sn  x  S  x  S  S"  is  a  state  transition  function  for  the  differential 
equation. 

We  now  consider  that  the  circuit  to  be  simulated  has  as  an  input  a  clock  with 
known  period  T  that  is  much  smaller  than  the  simulation  interval.  In  addition, 
we  assume  the  sequence  formed  by  sampling  the  state  at  the  beginning  of  each 

clock  cycle,  x(0),z(T),x(2T) . z(mT),. ...  changes  slowly  as  a  function  of  m, 

the  clock  cycle  number.  For  such  a  circuit,  a  designer  is  usually  interested  in 
the  transient  behavior  of  an  envelope  of  the  solution,  which  we  define  as  the 
continuous  function  derived  by  interpolating  the  sequence  formed  by  sampling 
the  state  every  time  interval  T.  Note  that  our  use  of  “envelope”  is  not  the  com¬ 
mon  usage.  Here,  the  envelope  is  not  unique  given  z(<);  the  envelope  generated 
by  interpolating  the  sequence  z(0  +  r),  z(T  +  r),  z(2T  +  r), . . .  depends  on  t. 

A  “differential-like”  equation  can  be  written  for  the  elements  of  the  sequence 
x(0),z(T),x(2T ), . .  associated  with  one  envelope  of  the  solution  to  (1).  Ap¬ 
plying  (2),  the  elements  of  the  sequence  can  be  related  by 


(r((m  —  1)7’),  (m  —  \)T,mT). 


(3) 


The  relation  in  (3)  indicates  how  rapidly  the  initial  point  of  each  clock  cy¬ 
cle  changes  from  one  cycle  to  the  next,  and  in  that  sense  is  like  a  differential 
equation.  This  similarity  can  be  exploited  to  derive  methods  for  approximately 
solving  for  the  z(mT)'«.  For  example,  the  value  of  z((m  4-  l)T)  can  be  approx¬ 
imated  by 


x((m  +  l)T)  -  z((m  +  1)7)  «  (/  -  l)d(z(m7),  m7,  (m  +  1)7),  (4) 

which  is  loosely  analogous  to  solving  a  differential  equation  by  forward-Euler. 

To  compute  an  envelope  for  a  system  with  period  7  using  a  forward-Euler 
envelope-following  algorithm  with  a  fixed  cycle-step  /,  a  simple  repetitive  two- 
step  process  can  be  used.  Given  z(0),  the  first  step  is  to  calculate  z(7)  by 


3 


solving  (1)  over  the  interval  (0,T)  using  a  standard  discretisation  technique. 
The  second  step  is  to  set  x(lT)  =  *(T)  +  (/  -  1)[*(T)  -  *(0)J.  This  process  is 

repeated  to  compute  x(2/T),  x(3 IT) . Note  that  calculating  the  solution  over 

a  long  interval  only  requires  solving  the  differential  equation  every  l*  cycle. 

Although  simple  to  describe,  a  forward-Euler  based  envelope-following  method 
is  not  very  effective  for  solving  switching  circuits  because  maintaining  stability 
severely  limits  the  size  of  the  cycle-step  /,  just  as  with  the  standard  forward- 
Euler  algorithm.  A  more  stable  algorithm  is  to  approximate  the  value  of 
x((m  +  /)T)  by 

x((m  +  l)T)  -  x(mT)  mt I*(x((m  +  /  -  1)T),  (m  +  l-  1  )T,  (m  +  l)T),  (5) 

which  is  analogous  to  baekward-Euler  for  the  differential  case.  This  approach 
allows  for  larger  cycle-steps  than  the  forward-Euler  based  approach,  but  leads  to 
more  a  complicated  equation  to  compute  each  cycle-step.  To  see  this,  consider 
computing  x(/T)  given  x(0)  based  on  (5).  An  x((l  -  1)T)  must  be  found  such 
that  when  used  as  an  initial  condition  for  (1),  the  x(lT)  computed  with  standard 
discretization  techniques  satisfies  x(lT)  —  x(0)  =  l[x(lT)  —  x((l  —  1)T")]  This  is 
a  boundary  value  problem,  and  is  in  general  difficult  to  solve.  For  the  case 
of  switching  power  or  filter  circuits,  the  above  boundary  value  problem  can 
be  solved  efficiently  using  a  Newton  method,  and  this  is  presented  in  the  next 
section. 

3  Solution  by  Newton 

As  mentioned  in  the  previous  section,  each  cycle-step  of  the  backward-Euler 
envelope-following  method  requires  the  simultaneous  solution  of 

x((m  +  l)T)  -  x(mT)  =  J[x((m  +  l)T)  -  x((m  +  /  -  l)r)]  (6) 


and 

x((m  +  l)T)  —  x((m  + 1  —  1)T)  = 

<t>(x((m  +  /  -  l)r),(m  +  /  -  l)r,(m  +  l)T).  (7) 

for  x((m  +  /  -  1)T)  and  x((m  +  l)T),  where  where  x(mT)  is  presumed  known. 
Therefore,  (6)  and  (7)  represents  2n  equations  in  2n  unknowns. 

An  iterative  Newton’s  method  can  be  applied  to  solving  the  above  system. 
In  general,  the  Newton  method  applied  to  the  problem  of  finding  an  x  6  3? 
that  F(x)  =  0,  F  :  8"  —  yields  the  iteration  equation  Jf(x*)[xt+1  -  x‘j  = 
-F(xk),  where  k  is  the  Newton  iteration  count  and  Jf  €  Sfnxn  is  the  Jacobian 
of  F.  Reorganizing  (6)  and  (7)  into  the  form  to  apply  Newton’s  method  leads 
to 

(I  -  f)x((m  +  l)T)  +  /x((m  +  /  -  l)T)  -  x{mT) 
x((m  +  l)T)  -  x((m  +  /  -  l)T)  -  <t>() 
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=  F(x((m  +  l)T),  x((m  + 1  -  1  )T))  =  0 


(8) 


In  this  case,  Jp(z((m  -t-  /)T),x((m  +  /  -  1)7’))  is  given  by 


f  (1  -  /)/„  //„  1 

l  In  -/«-&*()] 


(9) 


where  J„  is  the  identity  matrix  of  size  n. 

The  most  time-consuming  computation  in  this  Newton  iteration  is  evaluating 
Jr  and  F,  which  involves  computing  the  state  transition  function,  £(x((m  +  /  — 
1  )T),  (m  +  /  —  l)T,(m  +  l)T),  and  its  derivative.  The  state  transition  function 
can  be  evaluated  by  numerically  integrating  ( I)  from  (m  +  /  -  1)T  to  (m  +  l)T 
given  x((m  +  /  -  I)T).  The  derivative  of  the  state  transition  function,  referred 
to  as  the  sensitivity  matrix,  represents  sensitivity  of  x((m  +  f)T)  to  changes  in 
x((m  +  /  —  1)7’),  and  can  be  computed  with  a  small  amount  of  additional  work 
during  the  numerical  integration. 

To  show  how  the  computation  of  the  state  transition  function  and  its  deriva¬ 
tive  fit  together,  consider  numerically  integrating  (1)  with  backward-Euler, 
which  we  chose  for  its  simplicity  and  because  it  appears  to  be  one  of  the  bet¬ 
ter  formulas  for  clocked  analog  circuits.  Given  some  initial  time  to  and  some 
initial  condition,  x(t0),  applying  backward-Euler  to  (1)  results  in  the  following 
algebraic  equation, 


g(x(l0  +  A).  x(<0))  =  ^(p(*(<o  +  /»))  -  P(z(to)))  +  /(*(<0  +  *))  *  0  (10) 

where  h  6  9?  is  the  timestep.  Note  we  have  dropped  explicitly  denoting  the 
dependence  of  p  and  /  on  the  input  u  for  simplicity. 

Equation  ( 10)  is  usually  solved  with  Newton-Raphson,  for  which  the  iteration 
equation  is 

J,(x<‘)(to  +  h))(x<*+1>(to  +  h)  -  x<*>(to  +  h))  = 

-9(x(k){t0  +  h),x<k\t0))  (11) 

where  Jt(x(t))  €  $Nxfi  is  the  Frechet  derivative  of  the  nonlinear  equation  in 
(10)  and  is  given  by 


J,(x(t)) 


dg(z(t),  )  _  1  dp(x(t))  df(x(t)) 

dx(t)  h  dx(t)  flx(f)  ‘ 


(12) 


Solving  (10)  yields  an  approximation  to  x(!o  +  h)  =  ^(x(f0Mo.fo  +  h). 
Implicitly  differentiating  (10)  for  x(fo  +  h)  with  respect  to  x(fo)  yields 


+  h)) 


dx(t0  +  h) 
dx(t0) 


1  dpjx(t0)) 
h  dx(t0) 


(13) 
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Given  a  z( to).  (10)  can  be  repeatedly  applied  to  find  *(to+T)  =  ♦(*(to).  to.<o+ 
T),  and  (14)  can  be  repeatedly  applied  to  find  the  sensitivity  matrix  d*(fo  + 
T)/dx(t0)  =  d<j>(x(t0),to,to  +  T)/dx(to)  [aprille72].  Note  that  J /  is  required  in 
both  (12)  and  (14),  and  thus  the  sensitivity  matrix  update  can  be  made  more 
efficient  by  factoring  J,  once  and  using  it  for  both  computations.  However,  the 
sensitivity  matrix  is  still  expensive  to  compute,  because  it  is  an  N  x  AT  full 
matrix. 


4  Implementation  and  Test  Results 

An  envelope-following  method  has  been  implemented  in  the  NUsvrit  [lcundert88] 
simulation  program.  The  program  is  written  in  “C ”,  and  runs  under  the  UNIX 
operating  system.  The  program  uses  a  trapazoidal-rule  based  envelope-following 
algorithm  in  which  the  cyclesteps  are  selected  based  on  local  truncation  error. 
The  boundary  value  problems  generated  at  each  cyclestep  are  solved  with  the 
Newton  method  described  above. 

Two  techniques  are  also  used  to  improve  the  efficiency  of  the  basic  algorithm. 
Although  switching  power  supplies  and  filters  are  not  linear  circuits,  the  state 
transition  function  over  one  cycle  is  a  nearly  affine  (linear  plus  a  constant).  This 
property  can  be  exploited  to  reduce  the  computation  by  only  computing  Jr  f or 
the  boundary  value  Newton  method  on  the  first  iteration  at  each  timestep.  This 
is  a  significant  savings,  as  then  the  sensitivity  matrix  need  only  be  computed 
once  per  timestep.  The  second  technique  that  reduces  the  computation  time 
is  based  on  the  fact  that  the  explicit  forward-Euler  based  envelope-following 
method  is  stable  for  cyclesteps  of  size  two  or  less.  This  implies  that  if  truncation 
error  considerations  determine  that  a  cyclestep  of  two  or  less  is  appropriate,  the 
boundary  value  problem  can  be  avoided,  and  the  step  can  be  computed  easily. 

Exactly  how  the  envelope-following  method  behaves  can  be  seen  by  examin¬ 
ing  figures  1  and  2  below.  Figure  1  is  the  solution  to  a  buck-derived  circuit(from 
[hsiao87])  computed  with  classical  direct  method,  where  the  exact  solution  is 
only  shown  between  1  and  2  milliseconds  because  of  point  limitations  in  the 
plotting  program,  and  Figure  2  is  the  result  produced  by  the  envelope-following 
method.  As  can  be  seen,  the  envelope- following  method  computes  many  fewer 
cycles,  but  the  ones  computed  match  with  the  direct  method. 

In  the  table  below  we  present  a  comparison  between  the  cpu  time  used  by 
classical  and  envelope-following  methods  in  simulating  the  start-up  transient 
from  three  types  of  switching  power  supplies,  a  push-pull  flyback  converter, 
fly,  a  resonant  converter  (from  [casey87]),  res,  and  a  buck-derived  circuit,  buck, 
and  the  step  response  from  a  switched-capacitor  low-pass  filter.  In  each  case, 
the  clocking  is  provided  by  a  user-defined  source.  As  can  be  seen  from  the 
table,  the  envelope-following  method  can  be  very  efficient,  particularly  when 
the  simulation  interval  is  long  compared  to  the  clock  cylce. 
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Circuit 

Nodes 

Inteval/Clock 

Classical 

Env.  Follower. 

scop 

13 

200 

601 

153 

buck 

8 

1000 

940 

97 

quasi 

7 

200 

144 

38 

fly 

32 

40 

274 

167 

Table  1:  CPU  Time  (in  seconds)  Comparisons  for  Classical  vs  Enve¬ 
lope-Following  simulation,  based  on  a  SUN4 


5  Conclusions  and  Acknowledgements 


In  this  paper  it  is  shown  that  an  envelope-following  approach  to  the  simulation  of 
switching  power  and  filter  circuits  can  provide  substantial  speed  improvements 
over  classical  simulation  methods.  Several  aspects  of  the  method  are  still  under 
investigation.  In  particular,  it  has  been  observed  that  most  of  the  entries  in 
the  sensitivity  matrix  remain  close  to  zero,  and  how  to  exploit  this  is  being 
considered.  In  addition,  the  effectiveness  of  the  envelope-following  is  somewhat 
dependent  on  where  the  cycle  boundaries  are  placed,  and  an  automatic  selection 
method  is  desirable. 
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